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f-H \ We present a next-to-next-to-leading order (NNLO) realization of a general quark 

Qh' mass scheme (S-ACOT-x) in deep inelastic scattering and explore the impact of 

O NNLO terms on heavy-quark structure functions F% L (x, Q). An amended QCD fac- 

torization theorem for DIS is discussed that validates the S-ACOT-x scheme to all 
orders in the QCD coupling strength. As a new feature, kinematical constraints on 
| collinear production of heavy quarks that are crucial near the heavy-quark threshold 

are included in the amended factorization theorem. An algorithmic procedure is out- 
lined for implementing this scheme at NNLO by using mass-dependent and massless 
results from literature. At NNLO, the S-ACOT-x scheme reduces scale dependence 
of heavy-quark DIS cross sections as compared to the fixed-flavor number scheme. 
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INTRODUCTION 



In a modern global QCD analysis of parton distribution functions (PDFs), several factors 
rN | are comparable in magnitude to next-to-next-to-leading order (NNLO) radiative contribu- 

tions in the QCD coupling strength a s . Among these factors, dependence of QCD cross 
sections on masses of heavy quarks, m c and m&, can be significant. Global fits are sensitive 
to two types of mass effects, kinematical suppression of production of c and b quarks near 
respective mass thresholds in deep inelastic scattering (DIS), and large radiative contribu- 
tions to collinear production of cc or bb pairs at large collider energy. The first effect - 
suppression of DIS charm production near the threshold - must be carefully estimated when 
obtaining PDF parametrizations in order to accurately predict key scattering rates at the 
Large Hadron Collider jlj. The second effect is tied to an observation that c and b quarks 
behave as practically massless and indistinguishable from other massless flavors in typical 
Tevatron and LHC observables. It is therefore natural to evaluate all fitted cross sections 
in a "general-mass" (GM) factorization scheme, which assumes that the number of (nearly) 
massless quark flavors varies with energy, and at the same time includes dependence on 
heavy-quark masses in relevant kinematical regions. 

In this paper, we study NNLO quark mass terms in the default GM scheme of CTEQ 
PDF analyses called "S-ACOT-x" . Since its inception in 1993 0, the ACOT scheme has 



2 



undergone evolution based on the work in [3|-|5j. The S-ACOT-x version of the ACOT scheme 
is employed successfully to compute heavy-quark cross sections in recent NLO CTEQ6.6, 
CT09, and CT10 global fits 

The S-ACOT-y scheme is motivated by the QCD factorization theorem for DIS with 
massive quarks |3|, which provides the scheme's organizational backbone and key methods. 
In Sec. [Til we demonstrate explicitly how to amend the QCD factorization theorem in order 
to validate the S-ACOT-x scheme to all orders of a s . We then apply this scheme at NNLO 
to neutral-current DIS production, which provides the bulk of the DIS data, and for which 
all components of the calculation are readily available. 1 



Compared to other heavy-quark schemes available at NNLO |12Hl8|. our implementation 



aims to achieve more explicit analogy to the computation of NNLO cross sections in the zero- 



mass (ZM) scheme |19H21|. As another distinction, the S-ACOT-x scheme quickly converges 
to the fixed-flavor number scheme near the heavy-quark threshold as a consequence of the 
amended factorization theorem, without requiring supplemental matching conditions that 
are present in other general-mass schemes. 

In Sec. [Ill the S-ACOT-x cross sections are presented in the form that is reminiscent 
of counterpart ZM cross sections, up to replacement of some massless components by their 
mass-dependent expressions available in literature. This representation is based on a few 
compact formulas that include the desirable features existing in other heavy-quark NNLO 
calculations and lead to stable predictions for all scattering energies accessible. In Sec. Illfl 
numerical predictions are illustrated on the example of NNLO charm production cross sec- 
tions. They show that inclusion of the NNLO terms reduces theoretical uncertainties com- 
pared to NLO. 

Recent studies P, [22 24 j show that, at NLO, the LHC electroweak cross sections depend 
considerably on the mass scheme and parametric input for the charm mass m c in the PDF 
analysis, even though the combined HERA-1 data set (25| itself has a small total uncertainty. 
Yet, the upcoming combination of HERA-1 heavy-quark cross sections is expected to improve 
constraints on m c . The S-ACOT-x implementation brings theoretical predictions up to 
matching accuracy by including the NNLO terms. 



II. S-ACOT-x SCHEME: THEORETICAL FRAMEWORK 

Consider neutral-current DIS at energy that is sufficient to produce Ni light flavors (such 
as / = u, d and s) and one heavy flavor h ("charm") with mass rrih- We postpone the 
straightforward extension to production of several heavy flavors until Sec. Ill Dl 

The GM scheme is designed so as to enable quick convergence of perturbative QCD series 
involving heavy quarks at any momentum transfer Q. Perturbative QCD cross sections in the 
GM scheme must converge reliably near the heavy-quark production threshold (Q 2 « Tn^), 
as well as far above it (Q 2 ^> ml), and smoothly interpolate between the limits. When Q 
is of order rrih, it is most natural to include all Feynman subgraphs with heavy-quark lines 
into the hard-scattering function (Wilson coefficient function). Such approach is called a 
"fixed-flavor number" (FFN) factorization scheme. 0(a 2 ) coefficient functions for massive 



quark DIS production in this scheme have been computed in [26j-28]. At this Q, the NNLO 



1 For charged-current DIS, only massless 0, [Ic]] and some massive 11 1 NNLO coefficient functions have 
been computed. 
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coefficient functions in the GM scheme with Ni + 1 flavors are expected to reduce to the 
FFN massive cross sections in the FFN scheme with N flavors. On the other hand, at high 
virtualities (Q 2 > m 2 h ), the^NNLO GM cross sections should be indistinguishable from the 



NNLO ZM cross sections |19H21|. In this limit, the heavy-quark contributions are dominated 
by asymptotic collinear contributions that are also known to 0(a 2 ) fl2l. l29l-31|. 

A realization of such scheme called "ACOT" was developed in Refs. |2j, |32| and proven 
for inclusive DIS to all orders in Ref. A non-zero PDF is assigned in this scheme 
to each quark flavor that can be physically produced at the given Q. The more recent 
S-ACOT-x variant of the ACOT scheme adds two beneficial features. First, coefficient 
functions derived from Feynman graphs with initial-state heavy quarks are simplified by 
neglecting m c dependence 0, 0|. Second, threshold suppression is introduced by evaluating 
these coefficient functions as a function of % = x (1 + 4m 2 /Q 2 ) instead of Bjorken x 
Both modifications follow from the factorization theorem for inclusive DIS j3| and produce 
predictions that are simpler, yet numerically accurate. They are included as a part of the 
NNLO implementation that is now presented. 



A. Overview of QCD factorization 

A DIS structure function F(x, Q), such as F 2 or F L , is written in a factorized form as 

N f S N f i 

F(x, Q) = £ e 'E / | C ia ( x - 5, T± t as{fl) ) $ a/p(e , ^ 

i=l a=0 Jx * \^ V V / 

N}° Nf 

i=l a=0 

where $ a / p (£, /i) is a parton distribution function (PDF) for a parton type a, light-cone 
momentum fraction £, and factorization scale \x. Ci^ a (x,Q/n,mh/fi,a s (fi)) are Wilson coef- 
ficient functions evaluated at x = x/£. Convolution integrals over £ are indicated by 
Two sums appear on the right-hand side of Eq. flTJ, over all quark flavors % = l,...,Nj s 
that couple to the virtual photon with fractional electric charges = 2/3 or —1/3, and 
over parton flavors a in the PDF <3? a / p . The index a runs over quark flavors (a = 1, ...Nf 
for u,d,s,...) and the gluon (a = 0). Perturbative coefficients of neutral-current DIS are 
the same for quarks and antiquarks up to NNLO. For each combination of flavors i and a, 
summation of quark and antiquark contributions of these flavors is always implied, but not 
shown for brevity. 

Eq. ([T]) distinguishes between Nj s , the number of quark flavors produced in the final 
state (fs), and Nf, the number of active quark flavors in a s and PDFs. The distinction is 
important for the ensuing discussion, as generally Nj s is different from Nf [lj. Nj s is equal 
to the number of final-state flavors that can be produced at the given 7*p center-of-mass 
energy W = Qy/ljx — 1. All produced quark states can couple to the photon, so that the 
outer summation in Eq. (CQ) runs up to i = Nj s . 

On the other hand, Nf is a parameter of the renormalization and factorization schemes. 
It is commonly set equal to the number of quark flavors with masses lighter than Q. Only 
flavors with a < Nf have non-zero PDFs in the inner summation, but their actual number 
depends on the factorization scheme. 
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To determine Cj a , we calculate auxiliary structure functions for scattering on an initial- 



fs 



state parton b, F(e + b — > e + X) = Y^ii=i e l-^i,b- The coefficient functions Ci A are infrared- 
safe parts of Fif,- They enter convolutions together with parton-level PDFs $ a /b(^,yU,) for 
splittings b — > a, as 



Nf 



(2) 



0=1 



In the MS scheme, the parton-level PDFs are given by matrix elements of bilocal field 
operators that can be computed in perturbation theory. For example, the PDF for finding 
a quark q in a massless parton b, in the light-like gauge, is 



<M0 



dy 

2ir 



-£,p + y 



(6(p)|^(0,y-,0 T ) 7 + ^(0)|Kp)), 



(3) 



where the light-cone momentum components of the partons b and q are p^ = |p + ,0,0<r| 

and k^ = m 2 q /{2^p + ), 0t\, respectively, and p^ = (p° ±p 3 )/y/2. 

The functions Cj >a , and $ a /& can be expanded as a series of a s = a s (fi, Nf)/ (4tt): 



FifiW = FfllW + F$(x) + ai F$(x) + ..., 
C^(x) + a s CQ(x) + alC^(x) + ..., 

5 ab 5(l-0 + a s A^(0 + a 2 s A%(0 + - 



1(2), 



(4) 
(5) 
(6) 



In the last equation, (fc = 0, 1, 2, . . . ) are perturbative coefficients composed of DGLAP 

splitting functions Pf£\ such as = 2 -P^(£) In (/i 2 /m^) for the g hh splitting. 

By equating coefficients on both sides of Eq. (|2j), order by order in a s , we obtain 



F$(x), 

F$(x)- 

FS(x)- 



C 



(0) 



.4 



(i) 

ab 



UE - 



(i) 



.4 



(i) 

u b 



(7) 



Perturbative terms C\ k } in the coefficient functions are thus derived from perturbative ex- 
pansions for Fi^ and $ a /&, upon implied summation over the repeating index a. 

a (k) (k) 

The coefficients A ab in $> a /b consist of large or singular terms arising in F^ b when the 

(k) (k) 

momenta of a and b are collinear. Subtraction of convolutions of the A ab terms from F^ b 
on the right-hand side of Eqs. (J7]) produces a finite (infrared-safe) result. 

Depending on masses of a and 6, two forms of A^ b in these equations are possible. If 

(k) 

both a and b are massless, A ab contains a singular part, given in n = 4 — e dimensions 

by X^p=o(V e ) Ps p,a&> where s P)a b contains a DGLAP splitting function; and a finite part 

(logs+finite terms) of the form Ylp=o ^i^ 2 1 ^ 2 r) s 'p obi w here jJL is the factorization scale, and 
Hir is the parameter of dimensional regularization in the infrared limit. 
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When these "mass singularities" are subtracted as in Eq. ([7]), one obtains infrared-safe 
parts iy^ of iy* , denoted by a caret: 

The difference iy^ is finite, even though both the "bare" functions and the PDF coef- 
ficients A^j contain the singular l/e p terms, where p is a positive integer. 

If a massive parton a is produced from a massless parton b (as in g —> cc or u — >■ g — > cc) , 
the coefficients consist solely of logarithms and finite terms involving mass m a , 

A$UA) =X>V/™^UO- (9) 

The coefficient is finite for m a ^ 0, but — >■ oo when /i 2 3> m 2 ,. For such mas- 
sive quarks, the coefficients appear as subtractions from massive in expressions 
for C\ k ^. In accordance with the S-ACOT scheme, mass-dependent A^J only appear in 
explicit heavy-particle production, i.e., when transitions b(mt) = ())—>■ a(m a ^ 0) are in- 
volved. All other subprocesses use massless expressions, constructed from renormalized parts 

FiM given by Eq. ©. 

B. Heavy-quark component of inclusive F(x,Q) 

To construct the Wilson coefficient functions explicitly, we decompose F(x, Q) according 
to the (anti-) quark couplings to the photon [3]. Terms in which the photon couples to the 
light (/) or heavy (h) quark are designated as F\ and Fh, respectively: 

AT; 

F = Y J F l + F h , (10) 
i=i 

with 

F t (x, Q)=e 2 l Y J ® (?> Q)> (11) 

a 

and 

F h (z, Q) = c\ [Ch, a ® $ a /p] (s, Q). (12) 

a 

Note that this separation is purely theoretical: F\ and i 7 ^ cannot be measured separately 
or distinguished in another way. Furthermore, the heavy-quark component Fh is not the 
same as the semi-inclusive heavy-quark structure function Fh t si measured in experiments. 
The relation between F^si and Fh is clarified in Sec. Ill D| with the explicit formula given 
by Eq. flUD. 
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FIG. 1: Representative scattering contributions to F^x, Q). 



Concentrating first on the contribution with the photon coupled to h, we obtain its 
Wilson coefficients Ch, a {%) from the parton-level functions Fh ^ via Eqs. ([7]): 



(i) 



(2) 
h,h 



? (1) 
7 (2) 



hi 



c 
c 



(0) 

(0) 
h,h 

G 



h,h 



hg i h,l 
PS,(2) 



0; a 



(1) 



r h,h ^h,h 



A 



(i). 

hh i 



4(2) 
A /i 5 



.4 



(1) 
hg 



p (2) _ ^(0) 
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.4 

(2) 
hh 



hi 



U h,g 



A 



a 



(i) 

h,h 



A 



(1) 
hh 



(1). 



(13) 



(2) 

In these expressions, the coefficient / with the initial-state light quark depends on flavor- 



non-diagonal, or pure-singlet (PS), components of F^ and A^' , defined by 



I (2) 



G 



Gff + 5{jG 



I.J 



NS 



for G 



G 



(2) c,(2) 



'J 



and A 



(2) 



(14) 



(2) 

On the other hand, the coefficient C h h with the initial-state heavy quark depends both on 
the pure singlet (PS) and non-singlet (NS) components, as will be shown below. 

Representative diagrams for heavy-quark contributions in Eqs. (1T31) are shown in Fig. [U 
The reader may consult this figure frequently to identify various terms in the ensuing discus- 
sion. Propagators and external legs for quarks that are indicated by thick lines (thin lines) 
will eventually be computed with full mass dependence (in the massless approximation) . 
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The heavy-quark diagrams fall into two categories, those that do not involve a collinear 
approximation for scattering of heavy quarks (often called "flavor-creation", or FC, terms), 
and those that do (flavor-excitation, or FE, terms). While FC contributions must be evalu- 
ated exactly, the approximate nature of the FE terms allows some useful simplifications. 

At (9(a 2 ), the flavor-creation contributions include the coefficients Ff^, F^, and F^'^ . 
The heavy quarks appear in these terms inside the 0(a s ) Feynman subgraph for 7*^ — >• 
hh, connected by a gluon propagator to an initial-state gluon or a light quark. These 
contributions are evaluated with the exact kinematical dependence on m^, and hence are 
defined unambiguously. 

The FE cross sections are proportional to the heavy-quark PDF that approximates 
collinear production of heavy-quark pairs from light partons in the high-energy limit. Struc- 
ture functions and coefficient functions with an initial-state heavy quark, such as F^ and 

Cf^l, fall into this class. The FE coefficient functions reduce to unique MS expressions 
when m/j is negligible (3J, but, near the threshold, they may differ by non-unique powerlike 
contributions (mf l /Q 2 ) p with p > 0. Within the ACOT scheme, several conventions have 
been proposed to include the powerlike contributions in a way compatible with the QCD 
factorization theorem. 2 

Among these conventions, the "full ACOT scheme" Q evaluates the FE coefficient func- 
tions (Cffl, etc.) with their complete mass dependence. The simplified ACOT (S-ACOT) 

scheme 0, H| neglects all mass terms in and thereupon reduces tedious computations 
typical for the full ACOT scheme. The S-ACOT-% scheme Q adopted in our computation 
includes the most important mass dependence in and uses simpler zero-mass expressions 
everywhere else. 

If we denote the mass- dependent and massless quantities by uppercase and lowercase 
letters, and renormalized ZM functions by a caret, the S-ACOT-x convention for functions 
with initial-state heavy quarks is summarized as 

and 

^(f,f,f)=«(f,f,™, = 0)^ <?<!), (16) 

where 

*=»(i + $). m 

This form reflects an observation that the largest powerlike terms are associated with the 
constraint on the integration range in the convolution that reflects energy conservation in 
production of massive heavy-quark pairs (33|. This is practically achieved by evaluating 
ZM functions with \ instead of Bjorken x as an input parameter, and by constraining the 
momentum fraction in the FE channels to be in the range x < £ < 1- Justification for the 
X convention in the context of the proof of the QCD factorization is presented in Sec. Ill E[ 
Its numerical impact is discussed in Sec. IIII E[ 



2 The differences between these conventions are formally of a higher order in a s , but some conventions lead 
to faster perturbative convergence. 
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When the S-ACOT-% scheme is adopted, Eqs. (TT3j) become 

(18) 
(19) 
(20) 
(21) 

Lowercase functions on the right-hand side of these equations are given by ZM expressions. 

(1) ^(2) ^PS (2) 

Among all terms, only the structure functions F^ ' , F^ ■> and F hl are calculated with the 

exact mass dependence. The carets above F^ g and i*/ff indicate that the massless pole 

terms, Cg <8> Agg and <g) A g i , are subtracted from them. 

In the rest of the terms, the input longitudinal variable is set to be £ = x/C- The 
convolution of any such term f(x) with the PDF $(£) is 



c {0) 

c h,a 


= 6 ha 5(l 


- 2); 


c (1) - 


J h,hi 


rW 


- F w _ 

~~ h,g 


- 4 {1) - 

hgi 






^(2) 


- p(2) _ 
~~ h,g 


A (2) 


«& » 


A hgi 


(2) 
C h,h 


— J h,hi 


^(2) 


r h,l 


A PS, (2) 
~ A hl 



[f ® *] (0 



/'f (22) 



where C — X- [The naive massless approximation is ( = x.] 

The one- loop expressions /^2> F^, and A^ g can be found in jl, [H, . F^ g and F^f'^ 
coincide with the massive structure functions with initial-state gluons and pure-singlet light 



quarks in |26l . |27| . They are independent of A^. The expressions for A^ and A$ are 



computed as A^ and Ag q in [12| 

The 0(a 2 s ) contribution = f^ h corresponds to radiation of up to Ni + 1 flavors of qq 
pairs off an incoming quark h. It can be found as a sum of the pure-singlet and non-singlet 
ZM coefficient functions from Refs. (l9l-[2ll 35 . HI: 



(2) _ PS, {2) NS,(2) _ 7{2) _ -fPS,(2) 7NS,(2) / Q o\ 
c h,h ~ c h,h c h,h ~ J h,h ~ J h,h i" J h,h • 

In this equation, both pure-singlet and non-singlet parts of h are taken to be massless, 
which is one of the choices possible within the S-ACOT scheme. It is equally acceptable 

PS (2) 

to evaluate the pure-sing let F L h with full mass dependence, so that the corresponding 

coefficient function is given by the massive C$ in Eq. (J2Ij). In this case, F^'^ can be 
combined with the pure-singlet contribution with initial-state light quarks, also given by 

r {2) 
U h,l ■ 

The complete 0(a 2 s ) part of Fh(x, Q) then takes the form 

= e 2 h ® (<& fc/p + H, P ) + C*J ® E + Cg ® $ 9/p } , (24) 

where S(a;, //) is the singlet quark PDF summed over A^ = A^ + 1 flavors: 

E(z, fi) = J2 [®i/p( x > AO + A*)] • ( 25 ) 

We use Eq. (124")) for practical implementation of i^, with the coefficient functions computed 
as in Eqs. (|20|) and (121]) . 




(a) 



(b) 



(c) 



FIG. 2: Disconnected heavy-quark contributions to F^ S '^ (a,b,c) and A^j 3 '^' (d). 



(d) 

,NS,(2) 



C. Light-quark component of F(x,Q) 

By a similar argument, Eq. ([7]) serves as a starting point for finding coefficient functions 
for the light-quark component F[(x,Q). The corresponding Wilson coefficients are 
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») /I ( 2 ) 
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(i) 
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(1) 
/./( 

.(1) 



(1) 



o 



(1) 

,9 



.4 



(i). 

9' ' 



(26) 



The quark-to-quark Wilson coefficients Cj^f , G\^ and G\^, (for /' 7^ /) are decomposed into 

their PS and NS components as in Eq. (|14p . Non-singlet contributions F^ S '^ and A^ S 
contain squared matrix elements with heavy-quark lines that are disconnected from the 
initial-state proton, as in Fig. |2j These diagrams must be evaluated with full dependence on 
m/j. The rest of the coefficients in Eqs. (126]) do not contain disconnected heavy-quark lines. 
They are evaluated by ZM formulas. 

Explicitly, the non-singlet functions with mass dependence consist of two parts, arising 
either from Feynman diagrams with light partons only (designated as gu g ht), or with a heavy 
quark in the final-state emission or virtual loop (denoted by Gheavyi^h))'- 



(2) 



G — giight + G heavyi^h) , 

where G = F^'^ or A^f'^. The function gu g ht is evaluated in the ZM approximation, 
while G heavy ( m h) retains complete m h dependence. Masses can be neglected in the rest of 
Eqs. ( l26|) . so we get 



La 



/nf(2) _ /~fiv 
°Z,Z - 
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5iJ(l - x) 
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II, heavy 1 



f(2) 
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(27) 
(28) 
(29) 
(30) 
(31) 
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The 1-loop coefficients cf^ are well-known from Refs. (37M39|. The 2-loop massless contri- 
butions in Eqs. ([30]) and fl3T|) can be derived from the published ZM results according to the 
following procedure. Using the decomposition Eq. (I14p for in the ZM scheme, 

c hJ =c PS + 5 t3 c NS + 2 \ (32) 

where c PS and c NS are independent of the quark flavors i or j given that the masses are 
neglected, we write 

F(x, Q) = E,, a ef [ci, a ® $ a/p ] = e? {E, (c ps + ^c NS ) <g> $ i/p + c g ® $ 9/p } 

= [c™ ® + £^ { [<* ® E] + JV) [ Cfl S *„J } • (33) 

Here 

The singlet PDF S(x, /i) is given by Eq. (125]) . and the non-singlet sum of (anti-)quark PDFs 
is 



N 

,NS 



( 1 \ 

i=1 V / / 



Eq. (133]) expresses Ffx, Q) in the same representation as Eq. (4.1) in the N 3 LO calculation 
of DIS cross sections (36|. Comparing Eqs. (130]) and (13T]) with ZM coefficient functions in 
Section 4 of that reference (which are indicated here by an asterisk " * "), we find that 



,PS,m _ „(2,*) 

I,ps 



C ,w = C 



/iV>, (34) 



*£> = 4 2 9 *7iV/, (35) 
and 

^SS = cg2(n, = iVi) ) (36) 

with / = 2 or L for F 2 (x, Q) and F L (x, Q), respectively. 
The non-singlet heavy-quark coefficient function, 

<2»>.<?VmD = «f (2 '(x,Q7-D) + + §■» (J) <£>(*). (37) 

is composed of contributions of several classes shown in Figs. [2(a)- (c). Diagrams with 
real emission of a heavy-quark pair (as in Fig. 12(a)) in Fi(x,Q) contribute a function 

Lfq {x, Q 2 / m 1) m Eqs. (A.l) and (A. 2) of Ref. jloj. This contribution is combined with 
the virtual two- loop diagrams, cf. Fig. E(b), to produce the first term on the right-hand side 

of Eq. ([37]), in which L^ (2) (x, <2 2 /m|) is regularized by the plus prescription at x — > 1. 
Contributions with a heavy-quark polarization graph inserted into a one-loop 7*g scattering 
diagram, of the kind shown in Fig. [2(c), produce the second term in Eq. ([37]) . where cfi is 
available from Refs. 137-39 



In this derivation, we do not explicitly compute the virtual loop contribution in Fig. [2(b), 
but deduce it from the Adler sum rule 53-0. The sum rule states that the sum of the real 
and virtual contributions to F^ S ^ vy {x , Q 2 / 'ml) satisfies 

f KS vy MK)dx = Q. (38) 
Jo 
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With this rule, it can be demonstrated that the virtual contribution amounts to imposing 
the plus prescription on L^^ 2 \x, Q 2 /m 2 ) as in Eq. (137|) . 

In the asymptotic limit Q 2 ^> ml, F^^j® for the inclusive Fi contains large terms 

proportional to \\i(Q 2 /m 2 h ). Those coincide with the 0(a 2 ) non-singlet part A^f^ of 
the light-quark PDF that includes radiation of a heavy-quark pair as shown in Fig. E[d). 
^flMavy ls com Puted as A^ q S ^ 2 \z,m 2 / fi 2 ) in Eq. (B.4) of Ref. |l2j], which we evaluate as a 
function of of z = x/£ m accord with the S-ACOT-x scheme. 

When A^fj^ is subtracted from F^ S h ^ vy as in Eq. ( 130|) . the difference is free of the 
collinear logs. After the difference is combined with the light-quark-only contributions 
fi*Uight> we obtain the full non-singlet coefficient function C[ 2 i' NS , which coincides in the 
limit ml/Q 2 — > with its zero-mass MS expression in Eq. (8) of Ref. (lo|. For the longitudi- 

N S (2,) 

nal function Fi, the heavy-quark subtraction A { { he J vy is zero. Putting everything together, 
we obtain the final expression for the NNLO light-quark component, 

F[ 2) = e 2 {C?™ ® + <%) + c ps ^ ® S + eg S % /p } , (39) 
where the Wilson coefficients are listed in Eqs. (I29I) - (I31I) and (j34p - (l37p . 



D. Several heavy flavors, semi-inclusive heavy quark production 

Structure functions Fi and Fh in Eqs. fl39l) and (124j) are all that is needed to compute 
inclusive F(x, Q). Our expressions can be also extended to include two or more heavy-quark 
flavors: 

1=1 h=N l+1 

where the sum runs over all quark flavors satisfying 4m| < W 2 (i.e., up to the total number 
Nj s of quark flavors which can be produced in pairs for the given energy W). 

A clarification is in order that the heavy-quark component Fh of inclusive F(x, Q) (de- 
fined as the part proportional to the heavy-quark electric charge e|) is not directly measur- 
able. Rather, experiments publish the semi-inclusive (SI) heavy-quark structure function 
Fh,si(x, Q) that is determined from the cross section with at least one registered heavy me- 
son. In the case of F 2 at NNLO, F^si(x, Q) with h = c essentially coincides with the charm 
structure function F^ that is commonly measured by HERA experiments. 

While F^si must be defined with care to obtain infrared-safe results at all Q [l3|, for 
a global fit it is sufficient to approximate F^si in the following way (l8j|. At moderate Q 
values accessible at HERA, we define it as 

F h , SI (x,Q) = F h {x,Q) +5>^f' (2) ® + <%). (41) 

1=1 



Here Fh(x,Q) is the component with the heavy quark struck by the photon, cf. Eqs. (fT2|) 

N S ("2) 

and Q23J. Lj* w is the non-singlet part of the light-quark component Fi(x, Q) that contains 
radiation of a hh pair in the final state, as shown in Fig. Et^a). This is the same function that 
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was discussed in the previous subsection, and it is computed explicitly in Ref. However, 
since the virtual diagram in Fig. Mb) does not contribute to F^gj, the plus prescription is 

not imposed on L^f'^ in this case. 

The Fh^iix, Q ) fu nction that is thus defined is numerically stable in comparisons to 
the existing data [18|. Our numerical analysis shows that the contribution associated with 

TV S (2) 

Lj provides between and 3% of the semi-inclusive charm cross sections at Q < 10 
GeV, which is insignificant compared to typical experimental errors. 



E. Factorization and \ convention 

In the remainder of this section, we show that the S-ACOT-x scheme is fully compatible 
with the QCD factorization theorem for DIS. 

To see why the x convention is needed, consider again the heavy-quark contribution to 
F(x, Q) on the proton, 

F h {x, Q) = e 2 h V / ^C Kl Up, m h , Q-) * l/p (£, (42) 

This expression is the same as in Eq. (fl2"j) . but it assumes a low factorization scale /x ~ 1 
GeV, at which only PDFs for light parton flavors (I = 0, ...A r /) are present. The convolution 
limits in Eq. (142jl are such that the light parton I entering the coefficient function carries 
enough energy to produce one hh pair or more in the final state. As [i increases, a coefficient 
function Ch,h with an initial-state heavy quark is introduced at the switching point from TV; 
to N[ + i active flavors. This function does not automatically vanish outside of the physical 
range x < £ < !■ If Ch,h is defined so as to contribute in a wider range £ m ; n < £ < 1 at the 
switching point, with x < £, m i n < Xi then the same wider range will be preserved at all \x 
above the switching point by the DGLAP evolution. 

If W is the center-of-mass energy of the photon scattering on a light parton a, 

W 2 = (p a + q) 2 = Q 2 (Z/x-l), (43) 

it must be larger than the hh pair mass for any hh production subprocess to occur: 

Am 2 h <W 2 <W 2 = Q 2 (l/x - 1). (44) 

According to this condition, the scattering probability is non-zero only if the momentum 
fraction £ is in the range 

X < C < 1, (45) 

where x — x (1 + 4ml/Q 2 ) > x. 

If collinear approximations for flavor-excitation (FE) and subtraction terms in Ch, a violate 
this fundamental requirement, large spurious contributions from the unphysical kinematical 
region cancel to each order of a s , but survive as higher-order logarithmic terms. They can 
be eliminated by a supplemental condition that correct integration limits are always to be 
preserved, as in Eq. (T45]) . 3 



3 Even in the Q 2 3> m 2 h limit, convolutions with FE terms and subtractions could be in principle extended 
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We will now show how to apply this condition at any order by including it into the QCD 
factorization theorem. For this purpose, we examine the projection operator Z that encap- 
sulates the main rules of each factorization scheme It applies a set of approximations to 
Feynman graphs with leading momentum regions in order to enable all-order factorization. 

Feynman graphs containing the leading DIS contributions are composed of two-particle 
irreducible subgraphs H and T, joined by one parton line on each side of the unitarity cut. 
Each leading graph H ■ T involves integration over the momentum of the intermediate 
parton and summation over its spin components, 

H - T = E ( 46 ) 

a=g,u,u,d,d,... spins 

Virtualities of all momenta are of order Q in the hard subgraph H a (q, k), and they are much 
smaller than Q in the target subgraph T a (k,p). q^ and p^ are the photon's and proton's 
4-momenta. The nearly massless proton moves in the +z direction in the Breit reference 
frame. 

The purpose of the Z operator is to approximate the leading-power (logarithmically 
divergent) part of H ■ T by a simpler expression, denoted by H ■ Z ■ T, and to recast H ■ T as 

/d A k f dH . — 
J2^y J J2txY E Ha ( q > ^ Za ( 1 ' k] ^ T a(k,p) + non-leading power term. 



spins 



(47) 



The leading-power approximation H ■ Z ■ T provides the bulk of H ■ T. The non-leading 
power part is suppressed by terms of order 



/ highest virtuality in T 
V lowest virtuality in H 



with r > 0. 



When it is recursively applied to all leading subgraphs, the Z projection generates the 
factorized expression for the structure function, F = J2 a [C a ® f a ] + 0(Aqcd/Q)- Q 

The Z operation simplifies integration over the momentum and summation of the spin of 
the intermediate parton, and it also simplifies the hard graph H. The momentum Z M of the 
parton a that enters H is replaced by a simpler momentum P, e.g., P = t^p 11 if a is massless. 
The Z operation discards suppressed terms in H, notably, masses of the light quarks, and 
specifies how to handle the heavy-quark masses in H. 

Variants of the ACOT scheme are distinguished by the form of the operator Zh that 
acts on the H subgraphs with an incoming heavy-quark line. The target parts T a , and the 
operators Z\ and Z g for the H subgraphs with initial light-quark and gluon lines are the 
same in all variants. The PDFs in T a are defined by operator matrix elements as in Eq. ([3]) 
and retain mass dependence. 

The Zh operator is of the form 

Z h (l, k-V) = - A (2tt) 4 S h (T) S t 5(l + - T + )5(F - F)5 2 (r T ), 



to include contributions from < £ < x. This would not violate QCD factorization order by order, but 
destabilize higher-order terms. This is avoided by an implicit assumption that the FE convolutions in the 
ZM limit are restricted to the physical range x < £ < 1. 
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Scheme 


V* in Z h 


Sh 


m/i in H h (q,l) 


£ range in H h ■ T h 


ACOT 




T-y+m h 
?P+ 


m h ^0 


i (i + >/i + ^j<^<i 


S-ACOT 


(V,0,0 T ) 


7" 


m h = 


x < c < i 


S-ACOT-x 


(^l+4mVQ^°'" r ) 


7~ 


m h = 


x (1 + 4m2/Q 2 ) < £ < 1 



TABLE I: Components of the projection operator Zh(l, k; I) in three versions of the ACOT scheme. 



where Sjj(l) and St = 7 + are projectors on the leading spin components in H and T, 
respectively. The incoming heavy quark in Hh has an approximate momentum P, where 
the light-cone components of P in the Breit frame are / ± = \ l° ± Z 3 j /y/2 and = 0. With 
this representation, the Hh ■ Zh • Th integral assumes the form of a convolution over £, 



H h ■ Z h -T h = / y tr 



H h (q,l) 



S H (l) 



dk (Pkx 



tr 



(48) 



where £ = l + /p + is the ratio of large "+" momentum components. 

Table ffl collects expressions for P and S H (0 in the ACOT, S-ACOT, and S-ACOT-x 
schemes. It also specifies integration ranges in the Hh ■ Th convolutions and indicates if rrih 
is set to zero in the Hh subgraphs. We re-emphasize that the three schemes differ only in 
the Wilson coefficients with initial-state heavy quarks, i.e., in the flavor-excitation channel. 

In Table HJ we observe that the schemes differ by moderate terms in Hh{q,l) that are 
proportional to powers of m\/Q 2 . The form of these terms is not specified by the factor- 
ization theorem. It is allowed to simplify the FE hard-scattering contributions by setting 
rrih = as in the S-ACOT scheme. In the full ACOT scheme, the lower limit of integra- 
tion in Hh ■ Zh ■ Th is set by the kinematics of scattering of a massive quark into a massive 
quark, which violates the momentum conservation condition of Eq. ( |42|) for pair production 
of massive quarks from light-quark scattering. In the S-ACOT scheme, one is tempted to 
set the integration range to x < £ < 1, which is also incompatible with momentum conser- 
vation. One cannot just restrict the integration range to \ < £ < 1, as this disallows the 
lowest-order FE contribution c h °\ (g> h that contributes at £ = x. 

A better way is provided by rescaling £ — > £ = p + — > p + = p + / k, which leaves H(q, I) 
invariant, as it does not change / M , x = Q 2 / (21 ■ q), or other kinematical variables in H(q, I). 
At the same time, rescaling changes the integration range for £ in the convolution. 

Choosing k = (1 + Aml/Q 2 ), we obtain the S-ACOT-x scheme that has all desirable 
features: 

1. The proof of QCD factorization for the S-ACOT scheme in (3| also applies to the S- 
ACOT-x scheme, since the S-ACOT and S-ACOT-x schemes have the same H(q,l). 
The Z operation of the S-ACOT-x scheme upholds all expected properties that are 
listed in Sec. 9C of H. 



2. The integration over £ proceeds over the physical range x < £ < 1 in all channels. It 
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f f f 

FIG. 3: Integrands of convolution integrals F^&tg and A^®g with and without the \ prescription, 
plotted as a function of the momentum fraction £. 
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includes all physically possible scattering channels, but excludes kinematically prohib- 
ited £ values. 

3. The S-ACOT-x coefficient functions in the flavor-excitation channels are given by 
ZM expressions evaluated at x — x/£- Kinematical prefactors in front of the coefficient 
functions that are independent of £ are not modified. 

4. The target subgraphs T a , corresponding to the PDFs, are given by universal operator 
matrix elements that are the same in all ACOT-like schemes. 

5. When Q is much larger than rrih, the S-ACOT-x scheme reduces to the zero-mass MS 
scheme, without additional finite renormalizations. 

6. When Q is of order m^, the S-ACOT-x scheme approaches the FFN scheme faster than 
the S-ACOT scheme. Better matching on the FFN scheme is a consequence of the Z 
operation that satisfies energy conservation. It does not rely on propositions beyond 
the factorization theorem with energy conservation, such as conditions for derivatives 
of F(x,Q) [Hj or damping factors (la |. 



An illustration for the x convention 

The advantages of x rescaling can be demonstrated on the example of the 0(a s ) ^*g 
contribution, consisting of the gluon-initiated box graph and corresponding subtraction, 
and shown by the second and third graphs on the upper row of Fig. [T] 

[CS« 9 ](*,Q) = jf f s«,<3)if» (|) - J' f <Kf.O)42 (0 ■ < 49 > 

The integrands of the convolution integrals on the right-hand side, g{^,Q) F^ g {xlC) an d 

g(£,Q) ^.^(C/O) where ( = x or \i are plotted in Fig. [3]^a-c) as a function of £. The 
computation follows the numerical setup described in the next section. The scale on the 
£ axis is logarithmic: the convolution integrals are proportional to areas under respective 
integrand curves. For definiteness, we choose x = 0.01 and Q 2 =10, 100, and 10000 GeV 2 , 
but the same features are observed for other x and Q values. 

In the charm creation contribution, g®F^ , the integrand (red solid line) vanishes outside 
of the physical range \ < £ < 1 , where X ~ 0.018, 0.0108, and 0.010008 for Q 2 =10, 100, 
and 10000 GeV 2 . On the other hand, the naive choice £ = x of the S-ACOT scheme allows 
the g(£,,Q) ^^(C/O integrand (green dashed curve) to contribute in the unphysical region 
x < C < X- Its spurious contribution is comparatively large at the smallest Q. It is not 
fully canceled by the counterpart FE term c^ h ® c (x, Q) in the first upper graph of Fig. [H 
leading to a bloated higher-order uncertainty. 

The S-ACOT-x integrand g{^,Q) ^^(x/O vanishes below £ = x ( c f- the blue dash- 
dotted line). It is numerically moderate at physical £ values, £ > X? an d its integral cancels 

well with c^ h ® c (x, Q), as will be further demonstrated in Sec. IIV A[ Note also that the 

difference between the two definitions for the g(£, Q) A^ g (C/0 integrand is small in most of 
the physical range x < £ < 1- 
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As the virtuality Q increases, the difference between \ and x progressively reduces, letting 
£ to vary in a larger interval. Finally in (c), for very large Q, the S-ACOT and S-ACOT-x 
subtractions become identical. ®g approximates well the collinear splitting contribution 

that drives much of the shape of F^(x/09(0- When A^ ® g is subtracted from F^ g (g) g 
as in Eq. (149]) . it produces a moderate negative 0(a s ) contribution, which is further reduced 
at NNLO. These cancellations are further examined in Sec. IIVBI 



III. NUMERICAL EXAMPLES 

In this section, we show representative plots from our validation tests for NNLO inclusive 
structure functions F 2 (x,Q) and Fl(x,Q) computed according to the S-ACOT-x method. 
We focus on partial contributions in which the photon strikes a charm quark, given by 
Fh(x, Q) in Eq. (1121) for h = c, and for structure functions F = F 2 or F^. These contributions 
are referred to as F2 C (x, Q) and Fl c (x, Q) in the figures. The same comparisons have been 
repeated for bottom-quark functions Fib and Fh,, as well as for the full inclusive functions 

iv !+ i F h and alternative values of Bjorken x. The results of other tests 



show similar patterns and can be viewed at |45[. 

NNLO heavy-quark coefficient functions are computed using a program available from 



27j. This program tabulates two-loop heavy-quark coefficient functions in a form that allows 
fast evaluation of convolution integrals in the Q range covered by the experimental data. 
The PDFs in all comparisons are obtained by using the Les Houches Accord toy 



parametrization |46|, |47j at the starting scale Q = m c = a/2 GeV. Other input param- 
eters are a s {Qo) = 0.36 and the pole mass m c . 4 The switching between 3 and 4 flavors 
happens at Q — m c . The a s and PDFs are evolved to higher Q values by the HOPPET 
computer code (48|. 5 



A. Q dependence 

Fig. H] examines Q dependence of charm structure functions F 2c (left panel) and Fl c (right 
panel). They are computed to order a 2 s in all schemes, referred to as "NNLO" by the counting 
convention adopted here. The upper insets in both panels show predictions at x = 10~ 2 in 
the S-ACOT-x scheme, FFN scheme with Nf = 3, and ZM scheme with Nf = 4. The lower 
insets show ratios of the FFN and ZM predictions to the S-ACOT-x prediction. 

The left panel shows that the S-ACOT-x theory prediction for F 2c (x, Q) (blue solid line) 
is numerically close to the FFN prediction (red short-dashed line) at Q w m c and to the ZM 
prediction (magenta long-dashed line) at Q > 10 GeV. 

Similarly, in the right panel, the S-ACOT-x prediction for the longitudinal function 
Flc(%, Q) coincides with the corresponding FFN prediction at Q m m c and approaches the 



4 Our program can alternatively read MS masses as the input. In this case, the MS masses are later 
converted into pole masses, because the operator matrix elements are published as functions of the 
pole mass. 

5 Bottom-quark contributions are omitted in this comparison. The charm PDF is zero at Q < Qo = m c , 
but acquires a non- negligible value immediately above Qo through an 0(a 2 s ) discontinuity existing at the 
switching point. 
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ZM prediction at Q > 30 GeV. [Fl c (x,Q) is sensitive to mass-dependent corrections to 
scattering off longitudinally polarized photons. Its matching on the ZM prediction happens 
at higher Q values than in F 2c .] The S-ACOT-x prediction interpolates between the FFN 
and ZM predictions at intermediate Q values, precisely as expected. 
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FIG. 4: Comparison of F 2c (x,Q) (left) and F Lc (x,Q) computed at 0(a 2 ) in the S-ACOT-x (solid), 
FFN Nf = 3 (short dashed), and ZM Nf = 4 (long dashed) schemes, shown as a function of Q at 
x = 10~ 2 . 



B. Dependence on the factorization scale 

NLO computations leave substantial uncertainty in DIS charm-quark contributions due 
to the choice of the renormalization/factorization scale and differences in FE terms in the 
threshold region. NNLO terms drastically reduce these uncertainties. Factorization scale 
dependence, and its reduction from NLO to NNLO, is illustrated by Fig. [5j Reduction in 
uncertainties in the modeling of kinematical threshold effects is discussed in Sec. IIIIEI 

In Fig. E^a)-(c), predictions for F 2c (x, Q) in the S-ACOT-x and FFN (Nf = 3) schemes 
are plotted versus Bjorken x at representative Q 2 values of 4, 10, and 100 GeV 2 . The 
F 2c (x, Q) values on the y axis are multiplied by 10 3 a/x to better visualize the accessible x 
region. Central predictions are computed for /i = \/Q 2 + rn 2 , the default scale in heavy- 
quark DIS cross sections in the CT10 global analysis [6j. The error bands are obtained by 
varying the scale in the range Q < [i < \/Q 2 + Am 2 . 

At Q = 2 GeV in Fig. EJa), the NNLO S-ACOT-x central prediction (black solid line 
inside a green band) is slightly above the NNLO FFN prediction (short-dashed line inside a 
magenta band) and has smaller scale uncertainty than FFN. At Q below 2 GeV (not shown), 
the NNLO S-ACOT-x and FFN predictions get even closer. In contrast, the NLO S-ACOT- 
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FIG. 5: Comparison of predictions for F2 C (x, Q) in the S-ACOT-x scheme and alternative theoretical 
approaches at NLO and NNLO. Central predictions are for fj, = \J Q 2 + rn%, and the error bands 
are for Q < \i < \JQ 2 + 4m^. 
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X prediction (a long-dashed line inside a blue band) underestimates the NNLO FFN result 
and has wider scale dependence. 

As Q increases to 10 GeV (Fig. |5jc)), S-ACOT-x predicts more event rate than the FFN 
scheme both at NLO and NNLO. Altogether, the Q dependence in these figures is fully 
compatible with the matching of the S-ACOT-x results on the FFN and ZM results in the 
limits Q 2 « m 2 and Q 2 ^> m 2 , respectively. 



C. Alternative factorization schemes 

Figs. 15(a)- (c) also show NNLO predictions in alternative GM schemes, indicated by scat- 
tered symbols: the modified Thorne- Roberts (TR') scheme fl~34l6| as implemented by the 
MSTW'08 PDF analysis |H]; and FONLL scheme C used by the NNPDF collaboration pi. 
Their values are computed in the 2009 Les Houches benchmark study of GM schemes [50| 
by assuming the same x rescaling as the S-ACOT-x scheme. 

The three schemes are seen to be in good overall agreement, apart from minor differences 
traced to subtle variations in the NNLO implementations that the schemes provide. 

At Q = 2 GeV, the NNLO S-ACOT-x prediction lies slightly above the FONLL-C pre- 
diction and below the MSTW prediction. At Q = 10 GeV, the NNLO S-ACOT-x prediction 
becomes closer to the MSTW prediction and is still above the FONLL-C result. These dif- 
ferences can be understood by noticing that the compared schemes may differ in subleading 
perturbative terms. For example, the FONLL-C scheme includes a threshold damping fac- 
tor to match on the 3-flavor result near the threshold [l8|. The S-ACOT-x scheme is not 
using the damping factor and is expectedly close to the FONLL-C scheme at Q — > m c , but 
not strictly identical. The TR'/MSTW prediction includes a constant higher-order term (of 
order 0(aD) to improve smoothness of switching from 3 to 4 flavors at Q = m c . Neither 
S-ACOT-x nor FONLL-C includes such term, which is why they may predict smaller F 2c 
values at low Q. 



D. NNLO vs. NLO predictions 

Improved stability of the NNLO prediction in Fig. 0(a) can be appreciated by comparing 
it to the counterpart NLO result shown in Fig. |5Jd). Here, we collect a variety of NLO 
F 2c (x, Q) values at Q = 2 GeV, obtained in FFN, S-ACOT, TR'/MSTW, FONLL-A, and 
FONLL-B schemes. The S-ACOT and MSTW predictions are shown with the x scaling as 
well as without it. 

The spread in NLO values of i 7 2 C (x, Q) observed in the figure is extensive, nominally 
suggesting a large uncertainty in the resulting NLO PDF sets. However, when included in the 
PDF fits, the most extreme predictions for F2 C (x, Q) in this figure are excluded by the fitted 
DIS data, which prefer the values that are about the same as the (relatively unambiguous) 
NNLO result. In CT10 fit, the scale fi is set equal to *JQ 2 + m 2 , which brings the NLO 
S-ACOT-x prediction in agreement with the measured cross sections. Thus, according to 
the past global fits, the NLO cross sections can be reconciled with the heavy-quark data, but 
at the expense of tuning of the scale parameter, for each value of m c and rescaling variable. 
The key benefit of the NNLO calculation for i<2 C (a;, Q) is to automatically achieve such a 
good agreement, nearly independently of the factorization scale. 
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E. Threshold effects 



LH PDFs 0=2 GeV S-ACOT |_H PDFs Q=10 GeV S-ACOT 




(a) (b) 

FIG. 6: Dependence on the factorization scale fj, and the rescaling variable C(^)- Left: Q = 2 GeV; 
right: Q = 10 GeV. 



In the above discussion, our NNLO structure functions are computed using the optimal 
rescaling variable ( = x i n the FE heavy-quark contributions. The rescaling variable sig- 
nificantly improves convergence near the threshold, by excluding contributions to the FE 
convolution integrals that are kinematically impossible. 

Dependence on the rescaling variable can be explored with the help of a generalization of 
the x variable proposed in Ref. j33| and designated by (. The generalized rescaling variable 
( is implicitly defined by 

X = l + C A -(4m2)/Q2' ( 5 °) 

where A is a real number. Various choices of positive A produce a family of GM schemes in 
which C takes continuous values between x (no rescaling) and x (full rescaling) . Specifically, 
A = produces £ = x of the S-ACOT- x scheme, and A > 1 produces ( ~ x that corresponds 
to the plain S-ACOT scheme without rescaling. Negative A values (not shown) strongly 
suppress the FE contributions by setting ( > X- 

Fig. |6] shows the bands of variations in the NLO and NNLO S-ACOT i<2 C values at 
Q = 2 GeV (left panel) and 10 GeV (right panel) when /i and A are varied in the ranges 
Q < /i < \JQ 2 + 4m^ and < A < 100, respectively. If A is naively varied in the full range, 
at smallest Q values one obtains a large excursion in NLO predictions (light blue band), 
which is considerably reduced when going to NNLO (green band). As Q increases above a 
few GeV, dependence on A diminishes practically to nil, as in the right panel of the figure for 
Q = 10 GeV. Together with Fig. Fig. indicates that the physically motivated rescaling 
variable is more important at low Q and NNLO than the factorization scale choice. 
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IV. CANCELLATIONS BETWEEN FEYNMAN GRAPHS 



A. Cancellations at low Q 



In order to match on the FFN and ZM predictions, certain classes of Feynman diagrams 
inside the S-ACOT-x NNLO coefficient functions must cancel in the respective low-Q and 
large-Q regions. We will show how these cancellations come about in the case of the charm- 
quark function F 2c (x,Q), but the pattern holds for the bottom quark and other structure 
functions with suitable modifications. 

The cancellations are revealed by plotting differences between various matrix elements 
and collinear subtractions discussed in Section III which are established by applying the 
factorization formula at the parton level. 

In the Q w m c region, all FE contributions in Eqs. ffT8]) - (T2~Tj) must cancel to a high degree, 



in order for F2 C {x, Q) to reduce to the FFN matrix elements F^ 1 ^ and F^}- In the threshold 
region, the evolved charm PDF is effectively of order 0(a s ), 



<2) 



lim c(x, Q) ps a s (Q) 



< ® 9 



(51) 



a FE contribution to F 2c (x, Q) containing a coefficient c^ h is effectively of order 0(a k s 
Keeping this in mind, at order 0(a s ) the virtual-photon-charm scattering diagram with 
in Eq. (|18j) cancels the gluon-initiated subtraction term with in Eq. (119)) . and only the 
7*g fusion diagram F^J in Eq. (|T9|) survives in the total F 2c {x 1 Q). In this case, the difference 



D 



(i) 

c(°) 



JP) 



(x,Q) 



So) 

^h,h 



A 



(1) 
hg 



fj 



(52) 



where c(x, Q) and g(x, Q) represent the charm and gluon PDFs, must be close to zero. 

As the next order is included, the cancellation present in D^ 0) must further improve. 
Two differences quantify the cancellations to this order: 



D®(x,Q) = D%Ux,Q)-a 



(i) 



c(°) 



c(°>' 



.(0) 
^h,h 



A hg 



9 



{x, Q)-a 



So) 

^h,h 



A 



PS,{2) 
hi 



(x,Q), (53) 



in which the convolutions of cf^ h with 0{a 2 s ) operator matrix elements A^ and A^' (A> are 
subtracted from D^'^(x, Q); and 



(2) 



1^,(2) 



(x,Q) 



c { hl®^>®g 



(i) 



{x,Q), 



(54) 



which probes the cancellation between convolutions involving the coefficient c^ h . By plot- 



ting the difference D J 0) , we quantify how the 0(a s ) cancellation in D^L, proportional 



(i) 



to c 



(0) 
h,hi 



(2) 

improves upon inclusion of the NNLO corrections. The difference D c{1) quantifies 



convolutions with instead of c^\. 



an independent 0(a 2 s ) cancellation, of the same structure as D^ 0) , but composed of the 



)(i) n( 2 ) 

CC°)' c( ( 

provide visual guidance, these differences are compared to the FFN Nf 



The left panel of Fig. [7] shows x dependence of D^ 0) , D ( c(0 



and D^ {1) at Q 



2 GeV. To 
3 prediction at 
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0{a 2 s ) (solid black line), which is roughly equal to the total rate at this Q (cf. the previous 
subsection). We also plot the S-ACOT-x contribution of O(a s ) provided by c®\ (dashed 
blue line), nominally counted as the lowest-order contribution. While the LO contribution 
on its own is substantial comparatively to the FFN result, it is mostly canceled by the 
subtraction in Eq. (1521) . so that the resulting difference D^ Q) (long-dashed green line) is 
small. 

The cancellation in D^ 0) is further improved by including the next-order terms in D^ {Q) 
as in Eq. ( 1531) . The difference D c ' 0) (dot-dashed red line) and especially the counterpart 
difference D cm (dotted purple line) give decreasingly small contributions. They satisfy 



D 



(2) 



« 



D 



(2) 
C(°) 
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(i) 

C(0) 



<F 2>c {x,Q). 



(55) 



Therefore, as Q — > m c , the S-ACOT-% scheme exhibits an almost perfect match on the FFN 
computation by virtue of perturbative cancellations that improve with each order of a s . 



S-ACOT-chi, Q = 2 GeV 



S-ACOT-chi, Q = 10GeV 



FFNS nf=3 nnl 
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FIG. 7: Cancellations in the S-ACOT-x scheme at Q 2 « m 2 c (left) and Q 2 > m 2 (right). 



B. Cancellations at large Q 



A different cancellation pattern is observed when m c is negligible with respect to Q, when 
large logarithms collected in A^, etc. must be subtracted from massive contributions 
to obtain infrared-safe C^ k \ These cancellations are illustrated in the right panel of Fig. [7] 
by D^p and Dg 2 ^ quantifying collinear subtractions in contributions containing the 7*g box 
subgraph. The lowest-order difference Dg 1 ^ is equal to the convolution of the coefficient C^ g 



as defined in Eq. f|T9 



(x,Q) 



7 (i) 

h,g 



(x,Q) 



J0) 
~*h,h 



A hg 



(x,Q)\. (56) 
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In this expression the subtraction term matches on the 0{a s ) photon-gluon contribution 
represented by F^' . The x dependence of this matching is shown in the right panel of Fig. [7] 

for Q = 10 GeV. It can be seen that (blue short-dashed line) is quite small compared 
to the 0(a 2 s ) FFN result. 

The a?-order difference can be constructed as 



dW(x,Q) = D^\x,Q) + 



g V~>"W g 

which can be cast into the form 



{x,Q) + 



C 



(2) 
hi 



(x,Q) 



Df> = DW + al {F% « i + F^> « E - c£> » <» « , 
-4°i®4>S-ci°is< 5 - (2, ®E} (57) 

_____ '^('2') 

by virtue of Eqs. f[2"Uj) and (|2ip . At this order, the collinear logarithms arising in F^ ' g are 
canceled by c^' h ® and £g> A^ , and, similarly, the collinear term in F^'^ is removed 
by c i°i ® ^w' 5 ''' 2 '' • The ne ^ e ff ec ^ °f the subtractions is that D g 2) (the green long-dashed line) 
provides a small correction to D g l \ The perturbative series converge well for D g k ^: 

\Df>-D®\<\DV>\<F${x,Q). (58) 



C. Cancellations without kinematic rescaling 



Although the cancellations happen for any rescaling variable (, their perturbative con- 
vergence is slower for a non-optimal choice, such as ( = x. The differences D^ 0) , etc. for 
C = x are shown in Fig. [S] and, as one can see, they are generally larger than in the case of 
C — X- Nonetheless, the differences are reduced by going to NNLO, although not as fast as 
for the optimal rescaling choice. 



V. CONCLUSIONS 



We examined connections between multi-loop calculations for massive quark production 
and fundamental concepts behind QCD factorization. An NNLO calculation for neutral- 
current DIS with massive quarks was documented in a form that bears structural similarity to 



the NNLO computation in the zero- mass scheme |19|-|21|. This calculation is algorithmic and 



utilizes readily available NNLO expressions. The main formulas are presented by Eqs. (1401) , 
(EU), and 032). 

The conceptual foundation for the presented results is provided by the S-ACOT-x factor- 
ization scheme. The discussion emphasized several strong features of this scheme: its direct 
origin from the proof of QCD factorization for DIS relative simplicity, and compliance 
with phase space constraints on heavy-quark production at all energies. 

Throughout this study, we highlighted phenomenological importance of energy conser- 
vation in massive particle production. We have shown how the constraints from energy 
conservation can be satisfied in all channels as a part of the QCD factorization theorem. 
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S-ACOT (no rescaling), Q=2 GeV S-ACOT (no rescaling), Q = 10 GeV 




FIG. 8: Same as Fig. [?l in the S-ACOT scheme without rescaling. 



These constraints are included in the definition of the Z operation in the Collins' proof of 
QCD factorization by rescaling the partonic momentum fraction in flavor-excitation Wilson 
coefficients. The S-ACOT-x scheme thus realizes correct kinematical dependence solely by 
the means of the QCD factorization theorem and momentum conservation. 

Schemes of the ACOT family differ only in mass-dependent terms in heavy-quark Wilson 
coefficient functions. PDFs are given by the same operator matrix elements in all schemes, 
such as Eq. ([3]). Estimates of these PDFs from global fits converge to unique universal 
functions as order of the QCD coupling increases. Convergence is the fastest in the S- 
ACOT-x scheme. 

At NNLO, dependence of S-ACOT-x predictions on the factorization scale and other 
tunable parameters is reduced compared to NLO. Cancellations between classes of Feynman 
diagrams are stabilized once NNLO terms are included. The derivation of S-ACOT-x pre- 
dictions is simpler than in some other GM schemes, as it is carried out by assu ming a unique 



number of active flavors (A 7 /) and one set of universal PDFs in every Q range |12| . |29|. It is 
minimal, in the sense that it does not impose conditions on the Q derivatives of structure 
functions or introduce a damping factor jl8j]. Yet, after the NNLO terms are included, 
the S-ACOT-x predictions result in good agreement with the other GM schemes. As the 
default heavy-quark scheme of CTEQ PDF analyses, the S-ACOT-x scheme is going to play 
a crucial role in global fits at NNLO. 



Acknowledgments 

Many ideas in this paper were inspired by Wu-Ki Tung. We thank J. Smith for providing 
the computer code for computing massive NNLO DIS cross sections and for helpful comments 
on this work, and J. Huston, F. Olness, J. Pumplin, D. Stump, and other CTEQ members for 
useful discussions. MG and PMN appreciate stimulating communications with S. Alekhin, 
J. Bliimlein, A. Cooper-Sarkar, S. Forte, A. Mitov, J. Rojo, and R. Thorne. This work 
was supported in part by the U.S. DOE Early Career Research Award DE-SC0003870; by 



26 



the U.S. National Science Foundation under grant PHY-0855561; by the National Science 
Council of Taiwan under grants NSC-98-2112-M-133-002-MY3 and NSC-99-2918-I-133-001; 
and by Lightner-Sams Foundation. CPY appreciates hospitality of the National Center for 
Theoretical Sciences in Taiwan, where a part of this work was done. 



[1] W. K. Tung et al., JHEP 02, 053 (2007). 

[2] M. A. G. Aivazis, J. C. Collins, F. I. Olness, and W.-K. Tung, Phys. Rev. D50, 3102 (1994). 

[3] J. C. Collins, Phys. Rev. D58, 094002 (1998). 

[4] M. Kramer, 1, F. I. Olness, and D. E. Soper, Phys. Rev. D62, 096007 (2000). 

[5] W.-K. Tung, S. Kretzer, and C. Schmidt, J. Phys. G28, 983 (2002). 

[6] H.-L. Lai et al, Phys. Rev. D82, 074024 (2010). 

[7] P. M. Nadolsky et al, Phys. Rev. D78, 013004 (2008). 

[8] J. Pumplin et al., Phys. Rev. D80, 014019 (2009). 

[9] S. Moch and M. Rogal, Nucl. Phys. B782, 51 (2007). 

[10] S. Moch, J. A. M. Vermaseren, and A. Vogt, Nucl. Phys. B813, 220 (2009). 

[11] M. Buza and W. L. van Neerven, Nucl. Phys. B500, 301 (1997). 

[12] M. Buza, Y. Matiounine, J. Smith, and W. L. van Neerven, Eur. Phys. J. CI, 301 (1998). 

[13] A. Chuvakin, J. Smith, and W. L. van Neerven, Phys. Rev. D61, 096004 (2000). 

[14] R. S. Thorne and R. G. Roberts, Phys. Lett. B421, 303 (1998). 

[15] R. S. Thorne and R. G. Roberts, Phys. Rev. D57, 6871 (1998). 

[16] R. S. Thorne, Phys. Rev. D73, 054019 (2006). 

[17] S. Alekhin, J. Blumlein, S. Klein, and S. Moch, Phys. Rev. D81, 014032 (2010). 

[18] S. Forte, E. Laenen, P. Nason, and J. Rojo, Nucl. Phys. B834, 116 (2010). 

[19] J. Sanchez Guillen, J. Miramontes, M. Miramontes, G. Parente, and O. A. Sampayo, Nucl. 

Phys. B353, 337 (1991). 

[20] W. L. van Neerven and E. B. Zijlstra, Phys. Lett. B272, 127 (1991). 

[21] E. B. Zijlstra and W. L. van Neerven, Phys. Lett. B273, 476 (1991). 

[22] R. Thorne, PoS DIS2010, 053 (2010), arXiv:1006.5925. 

[23] S. Alekhin, S. Alioli, R. D. Ball, V. Bertone, J. Blumlein, et al. (2011), arXiv:1101.0536. 

[24] R. Placakyte (2010), talk at the PDF4LHC meeting, 

http:/ /indico.cern.ch/materialDisplay.py?contribId=:6&sessionId=2&materialId=slides&;confId= 103872. 

[25] F. Aaron et al. (HI and ZEUS Collaboration), JHEP 1001, 109 (2010). 

[26] E. Laenen, S. Riemersma, J. Smith, and W. L. van Neerven, Nucl. Phys. B392, 162 (1993). 

[27] S. Riemersma, J. Smith, and W. L. van Neerven, Phys. Lett. B347, 143 (1995). 

[28] B. Harris and J. Smith, Nucl.Phys. B452, 109 (1995). 

[29] M. Buza, Y. Matiounine, J. Smith, and W. van Neerven, Phys. Lett. B411, 211 (1997). 

[30] I. Bierenbaum, J. Blumlein, and S. Klein, Nucl.Phys. B780, 40 (2007). 

[31] I. Bierenbaum, J. Blumlein, and S. Klein, Phys. Lett. B672, 401 (2009). 

[32] M. A. G. Aivazis, F. I. Olness, and W.-K. Tung, Phys. Rev. D50, 3085 (1994). 

[33] P. M. Nadolsky and W.-K. Tung, Phys. Rev. D79, 113014 (2009). 

[34] E. Witten, Nucl. Phys. B104, 445 (1976). 

[35] S. Moch, J. A. M. Vermaseren, and A. Vogt, Phys. Lett. B606, 123 (2005). 

[36] J. Vermaseren, A. Vogt, and S. Moch, Nucl.Phys. B724, 3 (2005). 

[37] W. A. Bardeen, A. J. Buras, D. W. Duke, and T. Muta, Phys. Rev. D18, 3998 (1978). 



27 



[38] G. Altarelli, R. K. Ellis, and G. Martinelli, Nucl. Phys. B143, 521 (1978). 

[39] B. Humpert and W. L. van Neerven, Nucl. Phys. B184, 225 (1981). 

[40] M. Buza, Y. Matiounine, J. Smith, R. Migneron, and W. van Neerven, Nucl. Phys. B472, 611 
(1996). 

[41] S. L. Adler, Phys. Rev. 143, 1144 (1966). 

[42] G. Altarelli, Phys. Rept. 81, 1 (1982). 

[43] Y. L. Dokshitzer, G. Marchesini, and B. R. Webber, Nucl. Phys. B469, 93 (1996). 

[44] R. Brock et al. (CTEQ), Rev. Mod. Phys. 67, 157 (1995). 

[45] http://hep.pa.msu.edu/cteq/public/SACOTNNLO2011/. 

[46] W. Giele et al. (2002), hep-ph/0204316. 

[47] M. R. Whalley, D. Bourilkov, and R. C. Group (2005), hep-ph/0508110; 

http: / /hepforge. cedar, ac.uk/lhapdf/. 

[48] G. P. Salam and J. Rojo, Comput. Phys. Commun. 180, 120 (2009). 

[49] A. D. Martin, W. J. Stirling, R. S. Thorne, and G. Watt, Eur. Phys. J. C63, 189 (2009). 

[50] J. Andersen et al. (SM and NLO Multileg Working Group) (2010), p. 110, arXiv: 1003. 1241. 



